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^ , We review our proposal for a constructive definition of superstring, type IIB matrix 

■ model. The IIB matrix model is a manifestly covariant model for space-time and matter 

' which possesses A/" = 2 supersymmetry in ten dimensions. We refine our arguments to re- 

' produce string perturbation theory based on the loop equations. We emphasize that the 
space-time is dynamically determined from the eigenvalue distributions of the matrices. We 

• also explain how matter, gauge fields and gravitation appear as fluctuations around dynam- 

^ ' ically determined space-time. 
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O ■ §1. Introduction 

00 ■ 

o ■ 

Q"^ ' Several proposals have been made as constructive definitions of superstring the- 

. ory. ^) " The type-IIB matrix model, ^) " ^) a large N reduced model of maximally 

r~| I supersymmetric Yang-Mills theory, is one of those proposals. It is defined by the fol- 

"V* ' lowing action: 

^ I S = -^Tr (^^[A„A,][Af^,A-] + ^i;r'^[A^,i;]^ , (M) 

^ , where Afj_ and V are N x N Hermitian matrices, the former is a ten-dimensional 

^ ' vector and the latter is a ten-dimensional Major ana- Weyl spinor field respectively. 

It is formulated in a manifestly covariant way, which is suitable for studying nonper- 
turbative issues of superstring theory. Since it is a simple model of matrices in zero 
dimension, it does not possess degenerate vacua unlike its higher dimensional cousins. 
It is possible that the model possesses a unique vacuum, namely our space-time. If 
so, we can in principle predict the dimensionality of the space-time, low-energy gauge 
group and matter contents by solving this model. In such an endeavor, this model 
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can be studied by numerical simulations effectively. In this paper, we review the IIB 
matrix model and explain how the space-time appears dynamically and how the low 
energy gauge symmetry and difFcomorphism invariance emerge microscopically. 

We first list several important properties of the IIB matrix model. This model 
can be regarded as a large reduced model of ten-dimensional N = 1 supersym- 
metric SU{N) Yang-Mills theory. It was shown that a large N gauge theory 
can be equivalcntly described by its reduced model, namely a model defined on a 
single point. In this reduction procedure, a space-time translation is represented in 
the color SU (N) space, and the eigenvalues of the matrices are interpreted as the 
momenta of fields. Therefore, the basic assumption in this identification is that the 
eigenvalues are uniformly distributed. As a constructive definition of a superstring, 
on the other hand, we will see that wc need to interpret the eigenvalues of matrices 
as the coordinates of space-time points. The interpretation is T-dual to the above. 

Since our IIB matrix model is defined on a single point, the commutator of the 
supersymmetry which we inherit from the ten-dimensional J\f = 1 supersymmetric 
SU{N) Yang-Mills theory. 



(1-2) 



vanishes up to a field-dependent gauge transformation, and we can no longer interpret 

this supersymmetry as space-time supersymmetry in the original sense. However, 
after the reduction, we acquire an extra bosonic symmetry, 

5A^ = c^l, (1-3) 

whose transformation is proportional to the unit matrix 1 and an extra supersym- 
metry, 

fS('^A, =0 



(1-4) 



The linear combinations of these two supcrsymmetries (1-2) and (1-4), 

QW = QW + Q(2), Q(2) = _ Q(2))^ (1.5) 

satisfy the following commutation relations, 

[eig'(*) , e-2Q0)] = -2e-iJi,e2p''6^'^^ , (1-6) 

where is the generator of the translation (1-3) and i,j = 1,2. Therefore, if we 
interpret the eigenvahics of the matrices A^ as our space-time coordinates, the above 
symmetries can be regarded as ten-dimensional N = 2 space-time supersymmetry. 
Since the maximal space-time supersymmetry guarantees the existence of gravitons 
if the theory permits the masslcs spectrum, it supports our conjecture that the IIB 
matrix model is a constructive definition of superstring. This is one of the major 
reasons to interpret the eigenvalues of ^4^ as being the coordinates of the space-time 
which has emerged out of the matrices. 
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The second important and confusing property is that the model has the same 
action as the low-energy effective action of D-instantons. ^^"^ We should emphasize 
here the differences between these two theories, since we arc led to different inter- 
pretations of space-time. From an effective-theory point of view, the eigenvalues 
represent the coordinates of D-instantons in ten-dimensional bulk space-time, which 
we have assumed a priori from the beginning of constructing the effective action. 
On the other hand, from a constructive point of view, we cannot assume such a 
bulk space-time, in which matrices live. This is because not only fields, but also 
the space-time, should be dynamically generated as a result of the dynamics of the 
matrices. The space-time should be constructed only from the matrices. The most 
natural interpretation is that space-time consists of discretized points, and that 
the eigenvalues represent their space-time coordinates. Here we need to assume that 
the dynamics of the IIB matrix model is such that the resulting eigenvalue distribu- 
tions are smooth enough to be interpreted as Riemannian geometry. 

A final important property is that the typc-IIB matrix model has no free pa- 
rameters. The coupling constant g can always be absorbed by field redefinitions: 

/ ^ 9^/'^-^lx ('1.7) 

This is reminiscent of string theory where a shift of the string coupling constant 
is always absorbed to that of the dilaton vacuum expectation value (vev). In an 
analysis of the Schwinger-Dyson equation of the IIB matrix model, we introduced 
an infrared cut-off e, which gives a string coupling constant, ^st = l/Ne^. How- 
ever, through a more careful analysis of the dynamics of the eigenvalues, we have 
shown that there is no such infrared divergences associated with infinitely separated 
eigenvalues, and that the infrared cutoff e which we have introduced by hand can be 
determined dynamically in terms of A'^ and g. (The Schwinger-Dyson equation and 
the double scaling limit are discussed in §2.) 

We then explain several reasons we believe that the IIB matrix model is a con- 
structive definition of type-IIB superstring in addition to the symmetry argument. 
First, this action can be related to the Grccn-Schwarz action of a superstring by 
using the semiclassical correspondence in the large A'' limit: 

-^[, 

IV ^ y d^a^. (1-8) 
In fact, Eq. (1-1) is related to the Green-Schwarz action in the Schild gauge: 



(1-9) 



We need to integrate over the scale factor of the metric y/g in order to quantize the 
Schild action 



Z = J Dy/gDXDipe-^^'^^'''^. (1-10) 
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The matrix analog is the fohowing grand canonical ensemble: 



Z = J2 dAdi:e-^^^\ 



n=0 




If the large N limit is smooth, we expect that the (3 ^ (3c limit is identical to 
consider the microcanonical ensemble with fixed N and take N large. 

The correspondence can go farther beyond the above identification of the model 
with a matrix regularization of the first quantized superstring. Namely, we can 
describe an arbitrary number of interacting D-strings and anti-D-strings as blocks 
of matrices, each of which corresponds to the matrix regularization of a string. Off- 
diagonal blocks induce interactions between these strings. ^'^^ Thus, it must be clear 
that the IIB matrix model is definitely not the first quantized theory of a D-string, 
but a full second quantized theory. 

It has also been shown ^) that Wilson loops satisfy the string field equations of 
motion for type-IIB superstring in the light-cone gauge, which is a second evidence for 
the conjecture that the IIB matrix model is a constructive definition of superstring. 
We consider the following regularized Wilson loop:^)'^*^) 



Here, k!^ are the momentum densities distributed along a loop C; we have also 
introduced fermionic sources, A„. The symbol e is a short-distance cutoff of string 
world sheet. In the large N limit, e should go to so as to satisfy the double scaling 
limit. In Rcf. 2) it was proved, once we have taken the correct scaling limit, that 
the M = 2 supersymmetry is enough to reproduce the lightcone field equation of 
type IIB superstring from the IIB matrix model. In order to resolve the problem of 
the double scaling limit, we need to evaluate several quantities (i.e., an expectation 
value of the Wilson loop with almost zero total momentum) and there is still some 
subtlety as to how to take this double scaling limit in which we obtain an interacting 
string theory. It is discussed extensively in §2. 

Considered as a matrix regularization of the Green-Schwarz IIB superstring, 
the IIB matrix model describes interacting D-strings. On the other hand, in an 
analysis of the Wilson loops, the IIB matrix model describes joining and splitting 
interactions of fundamental IIB superstrings created by the Wilson loops. Prom 
these considerations, it is plausible to conclude that if we can take the correct double 
scaling limit, the IIB matrix model could become a constructive definition of type- 
IIB superstring. Furthermore, we believe that all string theories arc connected by 
duality transformations, and once we construct a nonperturbative definition of any 
one of them, we can describe the vacua of any other strings, particularly the true 
vacuum in which we live. 

The dynamics of eigenvalues, that is, dynamical generation of space-time was 
first discussed in Ref. 3). An effective action of eigenvalues can be obtained by 



r M 



w{C) = IV n exp{^e(A;^^A^ + A„V)} 



(M2) 



Ln=l 
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integrating all of the off-diagonal bosonic and fermionic components, and then the 
diagonal fermionic coordinates (which we call fermion zeromodes) . If we quench the 
bosonic diagonal components, {i = 1- ■ ■ N), and neglect the fermion zeromodes 

the effective action for coincides with that of supersymmetric Yang- Mills 
theory with maximal supersymmetry and vanishes respecting the stability of the 
supersymmetric moduli. The inclusion of fermion zeromodes as well as the non- 
planar contributions lifts the degeneracy, and we can obtain a nontrivial effective 
action for the space-time dynamics. In Rcf. 3) we estimated this effective action 
by perturbation at one loop, which is valid when all eigenvalues are far from one 
another, \xi — Xj\ :» yjg. Of course, this one-loop effective action is not sufficient to 
determine the full space-time structure, but we expect that it captures some of the 
essential points concerning the formation of space-time. One of important properties 
of the effective action is that, as a result of grassmannian integration of the fermion 
zeromodes, space-time points make a network connected locally by bond interactions. 
This feature becomes important when we extract diffeomorphism symmetry from our 
matrix model, which is discussed in §4. 

Once we are convinced that IIB matrix model is a constructive definition of su- 
perstring, we then have to give natural interpretation of low energy dynamics. That 
is, wc need to show how we can obtain local field theory in a low energy approxima- 
tion and the origin of local gauge symmetry in our space-time generated dynamically 
from matrices. We also have to show how the background metric is encoded in a low 
energy field theory in the space-time, especially the origin of diffeomorphism invari- 
ance. In Ref. 4) we have shown that, if we suppose that the eigenvalue distribution 
consists of small clusters of size n, the low-energy theory acquires SU{n) local space- 
time gauge symmetry. This gauge invariance assures the existence of a gauge field 
propagating in the space-time of distributed eigenvalues. Also, we have obtained 
a low energy effective action (a gauge-invariant kinetic action) for a fermion in the 
adjoint representation of SU(n), which becomes massless. The low-energy behavior 
for these fields is formulated as a lattice gauge theory on a dynamically generated 
random lattice, and hence supports our interpretation of space-time. We have also 
shown in Ref. 4) that the diffeomorphism invariance of our model originates in in- 
variance under permutations of the eigenvalues. Our model realizes the invariance in 
an interesting way by summing all possible graphs connecting the space-time points. 
The diffeomorphism invariance restricts the low-energy behavior of the model, and 
indicates the existence of a massless graviton in the low energy effective field theory 
on dynamically generated space-time. The background metric for propagating fields 
is shown to be encoded in the density correlation of the eigenvalues, while the dilaton 
vacuum expectation value is encoded in the eigenvalue density. A curved background 
can be described as a nontrivial distribution of eigenvalues whose density correlation 
behaves inhomogeneously. 

Both of these fundamental symmetries, local gauge symmetry and the diffeo- 
morphism symmetry, originates in the SU{N) invariance of the matrix model. It 
is quite interesting that these fundamental symmetries can arise from a very simple 
matrix model defined on a single point. These are discussed in §4. 
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The organization of this paper is as fohows. In §2, we review our analysis of 
the Schwinger-Dyson equation for the Wilson loops and discuss a problem on the 
double scaling limit. In §3, we briefly review our analysis on the dynamics of space- 
time, and show how a network picture of space-time arises. Here is an analogy with 
the dynamical triangulation approach to quantum gravity. We also discuss a recent 
result on numerical simulation. In §4, we discuss a possible origin of low energy gauge 
symmetry on a dynamically generated space-time and that of the diffeomorphism 
invariance. We also show that we can obtain a low energy effective action for several 
fields by using these low energy symmetries. Such low-energy effective theories are 
formulated as a lattice gauge theory on a dynamically generated random lattice. 
Section 5 is devoted to discussions. 



§2. Loop equations and scaling limit 



In this section, wc derive the light-cone string field theory of type IIB super- 
string from the Schwinger-Dyson equations for the Wilson loops (loop equations) . 
The purpose of this analysis is to verify that the IIB matrix model indeed reproduces 
the standard perturbation series of type IIB superstring and to fix how to take the 
scaling limit. Here we review the analysis in Ref. 2) with refinement. 

We regard the Wilson loop 

wiC) = Hiv{C)), 

viC) = Pexp (^ij^ \a{k^{a)A^ + X{(T)ij)j , (2-1) 

as the creation or annihilation operator for the momentum representation eigenstate 
of string lA;'*, X{(t)), where k'^{a) is a momentum density on the worldsheet and X{a) 
is its super-partner. We explain in §2.1 the reason why this interpretation is natural. 
The basic equations we consider in the following are the loop equations: 

= I dAd^-^{Tr{t-viC,))w{C2) ■ ■ ■ w{Ci)e-^}, (2-2) 

= J dAd^l^-^{Tr{eviCl))wiC2) ■ ■ ■ ^(Q)e-^}, (2-3) 

where is a generator of U {N) Lie algebra, and an equation which represents the 
local reparametrization invariance of the loop: 



2.L Wilson loops and light-cone setting 

In this subsection, we briefly sketch our basic idea for deriving the light-cone 
string field theory from the loop equations. For this purpose, let us consider only the 
bosonic parts. We emphasize here that we perform this simplification for explanation. 
In fact, as is explained later, we cannot obtain the light-cone string field theory from 
the bosonic reduced model. 
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We first explain a motivation to consider the Wilson loops like Eq. (2-1) in the 
following. Let us consider a gauge theory in a box with size a. We impose the 
periodic boundary conditions on the fields. Then a Wilson (Polyakov) loop is given 

by 



w 




where x'* is an arbitrary function which satisfies the condition 

xI'{2it) - x''(0) = ni'a. (2-6) 

Here the are the winding numbers in the )U-th directions. In the zero volume limit 
(a — 0), Eq. (2-5) reduces to 




The expectation values of the Wilson loops with nontrivial windings vanish if the 
translation {U{\)) symmetry is not spontaneously broken. This phase transition 
is known to be the deconfining transition in lattice gauge theory. This symmetry 
may be broken in the large N limit for bosonic reduced models. In order to obtain 
string theory in Minkowski space which is translation invariant, we have to keep 
the translation invariance by super symmetrizing the theory. We assume that the 
gauge theory is in the confining phase and hence well described by string theory. We 
further assume that there is no phase transition while we take a limit. Then the 
Wilson loop (2-7) must represent strings. Let us consider the Wilson loops with large 
winding numbers so that n'^a is finite in a — > limit. They represent the strings 
with no momentum but with many windings in the zero volume target space. In 
order to obtain the strings moving in the infinite-volume target space, we adopt the 
T-dual picture here. That is, we reinterpret x'^ as the momentum density, fc^, and 
obtain the (bosonic part of) expression (2-1). As is expected in the ordinary T-dual 
picture, the windings are converted to the total momenta, which is seen readily in 
the relation 

/■27r 

ni'a= / daki'ia). (2-8) 
io 

Thus we regard the Wilson loop (2-1) as the creation and annihilation operator for 
the momentum representation eigenstate of string {a) ^ \{a)) . Since is dual 
to the momentum in the expression (2-1), it can be interpreted as the space-time 
coordinate naturally. 

To make the connection with the light-cone string field theory, we consider the 
particular configurations of the Wilson loops, which we call the light-cone setting. 
The Schwinger-Dyson equations lead to the continuum loop equation as we explain 
shortly: 

(fc('')'^((7)2 + x'('')^((7)2)(^«(Ci) • • • w{Cb) ■ ■ ■ w{Ci)) = 0, (2-9) 
where for simplicity we consider only the free part and we denote i-gj^j^ as Xn{(7). 
We have also the local reparametrization invariance (the bosonic part of (2-4)), 

x'(^)^(c7)A;(f)(a)(^Ci)---^a)---MQ)) =0. (2-10) 
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We put k^^^^{a) = 1 for all the Wison loops by using the reparametrization in- 
variance so that we set the length of the strings to be equal to the + components 
of their total momenta. Wc consider the configurations of the Wilson loops which 
possess the identical light-cone time x"*". Namely we perform the functional Fourier 
transformations of the Wilson loop from k^^^~{a) to x^^^~^{a) and consider such con- 
figurations that x('')+(cr) = for all the Wilson loops. We also locate a group of 
the Wilson loops at = — oc which represent a particular initial state. After these 
prescriptions, we denote the Wilson loop by w{C): 

w{C) = w[x+{a) = x+,k+{a) = 1, k\ay,a = ~ p+] 

= J m-(a)e-^r^'^^^^('^)'="('^)^«(C)|,+(,)=i,,+ (,)=,+ , (2-11) 

where is the + component of the total momentum. This is the light-cone setting 
which is illustrated in Fig. 1. 

In the light-cone setting, the loop equation (2-9) reduces to 

k^'^-{a){w{Ci)---w{Cb)---wiCi)) 




(2-12) 



Fig. 1. The light-cone setting. We consider the configurations of the Wilson loops which possess 
the identical light-cone time a;"*". We also locate a group of the Wilson loops at — — oo which 
represent a particular initial state. We put fc+(a) = 1 for all the Wison loops by using the 
reparametrization invariance. 



IIB Matrix Model 



9 



Equation (2-10) reduces to 

x'W-(a)(ii)(Ci)---ii)(a)---^(Q)) 

= x'^''^\a)k^p)i{a){w{Ci) ■ ■ ■ w{Ck) ■ ■ ■ w{Ci)). (2-13) 

Using these equations wc can deform the Wilson loop locally from the constant x'^ 
surface and shift the value of k^'^^^{a) locally. Therefore we can recover general 
configurations of the Wilson loops. This fact ensures us to consider the light-cone 
setting. Integrating Eq. (2-12) over a and summing up over 6, we obtain an operator 
which shift the constant value of x+, that is, the light-cone Hamiltonian: 

i^{w{C,)---w{C,)---w{Ci)) 

= J2[ ^4(fc(')'(^)' + ^'(')n'^)')(^(Ci)---*(C6)---^(C0), (2-14) 
b •' ^ 

which is identical to the free light-cone Hamiltonian of bosonic string. Note that we 
can obtain a closed system of equations within the light-cone setting even though we 
consider the particular configurations of the Wilson loops. Our procedure is anal- 
ogous to obtaining a Hamiltonian in super-many-time theory, where a Hamiltonian 
density is defined on a general space-like surface. We consider a constant time surface 
and integrate a Hamiltonian density on it to obtain an ordinary Hamiltonian. 

2.2. Loop equations and loop space 

In this subsection, we represent our loop equations by the loop space variables 
and put them in the light-cone setting. 
The loop equation (2-2) is evaluated as 

= ^ ((Tr {[Au, + V}) v{Ci))w{C2) ■ ■ ■ ^(Q)) 

+z dafc«(a)(Tr(Pe^r'^'^'('='^'''(-')^''+^^^'('^'W) 

+i ^ dakP ia){Tr{Pe^C <i<^'(fc(')''(<^')A.+A(^)(.'W„(^^) 

6=2 •'^ 

^Pe'i: ci."(.(^V(."M,+A(^)(."W . . . ^^c,) . . . wiQ)), (2-15) 

where w{Cb) implies the absence of the Wilson loop w{Cb)- The first term in 
Eq. (2T5) represents the infinitesimal deformation of the string coming from the 
variation of the action S while the second and third terms represent the splitting 
and joining interaction of strings respectively. Similarly, the loop equation (2-3) is 
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evaluated as 



= ^{{Trir^'iA^, i;]v{Ci))w{C2) ■ ■ ■ w{Ci)) 



6=2 



6=2 

X 

We need to represent these equations as differential operators on the loop space. 
Then we treat field insertions in the loop such as 

Tt{[A^,^]v{C)). (2-17) 

We have the following identity, 

Tr{[A^,^]v{C)) 
( S 




- 77^ - 77^ TtL ( nO)^ + A(0)^ ) .,(0 



+higher order of e. (2" 18) 

Though the right-hand side of Eq. (2T8) vanishes in naive e ^ limit, we expect 
here that for finite A'' an ultraviolet cut-off e appears naturally on the worldsheet. 
We assume that e converges to zero in A/" — ^ oo limit such as e ~ (6 > 0). We 
should keep e and other quantities depending on A'^ in the process of the calculation 
and take ^ oo (e ^ 0) limit as the continuum limit in the final stage. 

We expect naively that we may ignore the terms except the first one in the 
right-hand sides of Eq. (2-18) since they are the higher order terms with respect to e. 
However, in general, they contribute to the renormalization of the lower-dimensional 
terms because in the loop equations they generate the divergences from operator 
product expansion. The ways we represent the loop equations are not unique since 
the right-hand sides of Eq. (2-18) can be expressed in infinitely many ways. However 
we expect that they are unique in the e — > limit in the loop equation. This is the 
universality of the differential operators, which is guaranteed by a power counting 
and symmetry as seen in §2.4. Here we may draw an analogy with the quantum field 
theory on the lattice. The lattice action may be expanded formally in terms of the 
lattice spacing a. Although the operators which arc suppressed by the powers of a 
formally vanish in the continuum, we cannot simply neglect them because they may 
renormalize the relevant operators. In fact we can write down many lattice actions 
which possess the identical continuum limit. 
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In the following, we often show only the naive leading terms in the loop equations. 
Note that we cannot, in fact, neglect infinitely many higher order terms of e as is 
discussed above. 

By multiplying Eq. (2-15) by -ik^^\{)), we obtain 

^ 111 U 6 S V 



g'^W e I (5A;(i)m(|) <5A;(i)m(-|) 



-I 

r2iT 



+ 







daA;(i)^(0)fc«(a)(Tr(Pe^r'^-'('='^^''('^')^''+^^^'(-')'^)) 
+ Y: r dak^'^^mP (a)(Tr (Pe^ ^ 'i'r'(ki''H<'')A,+m.'W 

(,_n Jo 



6=2 

XV 



(2-19) 



Similarly, Eq. (2-16) is represented as 



6=2 



6=2 

X 



(2-20) 



Here let us proceed to the light-cone setting. The + component of total momen- 
tum carried by the w[x~^{a') = , k'^{a') = 1, k^{(j'), \{u')] a' = ~ a)] is equal to 
a. Therefore we expect naively 

{w[x+{g') = x+, k+{a') = 1, k\a'), X{a)- a' = ~ a)\) = (2-21) 

because of the momentum conservation if o" 7^ 0. We also expect the same case for 
{w[x^{(j') = x"*", k'^ip-') = 1, A;*((7'), A((7'); a' = — a ^ P^\)- However these are not 
actually the cases since the eigenvalues of distribute in a finite range for finite 
iV, which violates the momentum conservation slightly. Therefore these two w have 
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support for small a generally. We define a nonvanishing quantity / by 

1= r da{w[x+ia')=x+,k+{a') = l,k\a'),X{a);a' = 0^a] 
Jo 

+w[x^{a') = x+, k+{a') = 1, F(a'), A(cj); a' = p+ - a ^ p+]). (2-22) 

We expect that / diverges as A'' — > oo and assume that its large behavior is 
I ~ N°- (a > 0). Then a part of the splitting interaction term in the loop equa- 
tions contribute to the first term in the loop equations representing the infinitesimal 
deformation of the string. 

Thus in the light-cone setting Eq. (2-19) leads to 

k^'^-{0){w{Ci)---w{Ci)) 



b=2 

X 



(2-23) 

where Z is defined by 

z = 2 (^1 + -^^e'MOaio^j (2-24) 

and Tr means that the corresponding Wilson loop is arranged in the light-cone 
setting. Here x^^, 9a and 6a are defined as follows: 

5 



Oa{<y) = J^y (2-25) 



Equation (2-20) is decomposed into two equations as follows. 
A«(0)(.S(Ci)---^(Q)) 
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^(V2a;'(^)-(0)^,(0) + x'^(i)(0)7L^«(0))(^S(Ci) • • • w{Ci)) 



^ 6=2 

p^./;<i."(fe(^)''K)A,+A(^)K)^)\^(^^) . . . . . . w{Q)), (2-26) 



X 

I rv(''^+ rJb) 

E 

6=2 

X ' 

= ^(x'W^(0)7L)-'|a«(0)(^«(Ci) • • • w{Ci)) 



xTt(Pe^/:'''^ d."(fca)''K)A,+A«(."W j^^^^) . . . ^^ci 

6=2 

^p,i/;<i."(fe(^)''K)A,+A(^)KW)^(C^) . . . ^(C7,) . . . w{Q))^. (2-27) 

We eliminate a half of fermionic degrees of freedom by using the above two equations 
just like eliminating the half of the fermionic degrees of freedom in the light-cone 
field theory by using the equation of motion. We can also rewrite Eq. (2-4) as 

x('^-{a){w{Ci)---w{Ci)) 

(2-28) 

Our task is summarized as follows. By using Eqs. (2-23), (2-26), (2-27) and 
(2-28) iteratively and repeatedly, we eliminate Xnd, dna x ~ , and evaluate the 
light-cone Hamiltonian, 

H = Y, j dak^^'^~{a). (2-29) 

In this process, various interaction terms of order l/I^ ^ 0{1/N''"') are generated, 
which represent processes where the k + 2 strings interact at one point, i.e., (k + 2)- 
Reggeon vertices. In this procedure, we should take the continuum limit by keeping 
the higher order terms of e vanishing in the naive e — ^ limit. This task is in general 
very hard to perform because we should treat infinite series. However we will discuss 
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in §2.4 that the continuum limit is completely controllable by an analysis based on a 
power counting of e and the symmetries. We should consider the supercharges first 
rather than the Hamiltonian in this more rigorous treatment. In the next subsection, 
we first sketch how we can derive light-cone Hamiltonian for type IIB superstring 
from the loop equations. 

2.3. Light-cone Hamiltonian 

Here we neglect the interaction terms in Eq. (2-23) and concentrate on the free 
part of light-cone Hamiltonian. We can see that if we ignore the non-quadratic terms 
and set £^ ~ a'^ ^ const, the free part of Eq. (2-23) has the same form as the free 
light-cone Hamiltonian of type HB superstring in the naive continuum limit except 
lacking the A^Aa term. In the following, we verify that we indeed obtain this term 
when we eliminate 9d in the free part of Eq. (2-23) using Eq. (2-27). Let us consider 
the naive leading contribution in the free part of Eq. (2-27): 

'W( \ i 2t 

e^^\a){wiC^) . . . w{Ci)) = %^JJ^3—Xi^){a){w{C^)w{C^) • • • w{&)). (2-30) 

We assume first that the free part of Eq. (2-23) correctly describes the free part of 
light-cone Hamiltonian. This assumption can be justified by showing that the non- 
quadratic terms are indeed negligible in the continuum limit except for finite renor- 
malization of the quadratic terms. Since we are dealing with free two-dimensional 
field theories, we can use standard techniques of conformal field theory to estimate 
the effects of the non-quadratic terms. We note that the xl^\a)'^ is of order l/e^ 

since we have a cutoff length e. So wc may expand x^^\aY = \/a€^+ : x^^^a)"^ :, 
where : y : denotes the normal ordered operator constructed out of y. The symbol 
a in the denominator is a quantity proportional to ^Jej g^l on dimensional grounds. 
In this way, the left-hand side of Eq. (2-30) becomes 

e'^^;^'^(^)7L(l - ■■ 4'\<^f ■■+■■ •)Ai^)(^)(^(Ci) • • • w{Ci)). (2-31) 

For example, a typical non-quadratic term in the free part of Eq. (2-23), 

6^fc(i)ncT)^:r'(i)^(a)7i„A?)(a):r'(^)na)7|,A;(^)(a), (2-32) 

generates the \a'\(j)\^a\a) term due to the order 1/e^ divergence from the operator 
product expansion of k^^'>'^{uYx'^^^\u)x'^^'>^{a). 

In this way, we see that non-quadratic higher-dimensional operators do not ap- 
pear in the continuum limit due to the suppression of the powers of e and only 
renormalize finitely the quadratic operators. Therefore we expect to obtain the fol- 
lowing free Hamiltonian in the continuum limit: 

(2-33) 
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Here we assume that terms with the negative powers of e and other finite terms 
such as fe'x" do not remain, which is guaranteed by an argument based on a power 
counting and symmetries as is seen in the next subsection. In the bosonic model, 
this is not guaranteed and the divergence remains in general. This is the reason 
why we cannot obtain the hght-cone string field theory from the bosonic reduced 
model. The Hamiltonian (2-33) is identical to that of type IIB superstring theory 
if Aa(o") and 6a{cr) are rescaled appropriately and rotated by a complex phase factor 
•q = exp(^) as follows: 



We elaborate more on this point in connection with the supercharges in the next 
subsection. Here we cannot determine the coefficients c, c', c" and c'", which include 
the effects of renormalizations, and in the next subsection we can fix them by using 
M=2 supersymmetry. Prom Eq. (2-33) we find a'^ ~ 9^1 and hence we should 
obtain the prescription of the scaling limit, g^N"'^^ ~ a!"^ ~ const. 

2.4. General proof 

In this subsection, we give a general proof of our assertion that the light-cone 

string field theory for type IIB superstring can be derived from the loop equations 
of the IIB matrix model. We use a power counting and a symmetry analysis based 
on M=2 supersymmetry, 50(8) invariance and the parity symmetry on the string 
worldsheet. 

2.4.1. Power counting and parity symmetry 

In order to perform a power counting for e, we first introduce a mass dimension 
on the worldsheet through the relation [e] = — 1 and determine the dimension of each 
field. The IIB matrix model action (1-1) is decomposed into 



By demanding the IIB matrix model action (2-35) to be dimensionless, we obtain 



K{ct) ^ r?A„(a), ^„((7) ^ r/*^„((7). 



(2-34) 




+ i(V2Va[^+, V'a] - V'aTidl^', V'd] ' ^alLi^' Aa] + V2^d[^-, V'd]) j • 

(2-35) 




(2-36) 




k-{a)A+ + k\a)A' - iXa{a)^Pa - iXa{<^)^a)) 



(2-37) 



Prom this, we also read off the relations 



[k^] + [A-] = 1, 
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[k-] + [^+] = 1, 

m + [A^] = 1, 

[A„] + [Va] = 1, 

[Xa] + [V'd] = 1. (2-38) 

Noting that we should set [k~^] to be zero since k~^{a) = 1 in our hght-cone setting, 
we can determine the dimensions of all quantities as follows: 

[k+] = 0, [k-] = 2, [k^] = 1, [A„] = i [A,] = I, 

[X+] = -1, [X-] = 1, [x'] = 0, = I, [0a] = 

[A+] = -1, [A-] = 1, [A^] = 0, [V.,] = ^, [i^a] = -\. (2-39) 

Next we define a symmetry which corresponds to the parity on the string world- 
sheet. It is seen easily that the IIB matrix model action (1-1) is formally invariant 
under the following transformation: 

ip -iip^. (2-40) 
This transformation flips the direction of the Wilson loop in the following way: 

w{C) = TV (pe'ft Mk^{'r)A,+Ha)i,)^ 

^ Tr JJ e^^o da{ki'{p+-a)A;j,-i\ip+-a)'tlj)^ _ (2-41) 

Therefore our theory has a symmetry under the transformation 

fc'^(cT) ^ kf'ip-^ - a), 

X{a) iX{p+ - a), (2-42) 

which wc identify with the worldshcet parity. We also obtain the parity transforma- 
tion for the dual variables and On- 

e{a) -ie{p+ - a). (2-43) 

2.4.2. M=2 supersymmetry 

We denote the supercharges generating the transformations 



[d'-^^ip = ^ri^''[Ai^,A'']e 

and 

fS^'^A, =0 



(2-44) 
(2-45) 
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by Q^^^ and Q^"^^ respectively. We can determine the dimensions and parities of the 
parameters e and ^ in Eqs. (2-44) and (2-45) by comparing their both sides, 

[Sa] = l^a] = ^, [£a] = [Cd] = ~, 

' ''S'aj £d ^ *£d) 

^ -iCa, ^ -i(d- (2-46) 

This fixes the dimensions and parities of the supercharges Q^^^ and Q^^^ since £aQa^^ + 
SaQi^"* +CaQa^'* + CdQ^^^ generates the transformations (2-44) and (2-45): 

m = m = -\, \(^'] = [Qf ] = \, (2-47) 

Q'i'^^Q^\QT^^QT■ (2-48) 

Here we note that Eqs. (2-47) and (2-48) are consistent with the anti-commutation 
relations 

{Q«,Q«} = 0, {Q(^),Q(^)} = 0, 

{g«,Qf} = V2F5,^. (2-49) 

2.4.3. Free parts of supercharges and Hamiltonian 

The supercharges Q*^^^ and Q*^^^ can be expressed as differential operators on the 
loop space using the Ward identities. In principle we can eliminate the operators 
k~ , X ~ , \a and 6a by repeatedly using the loop equations and the reparametrization 
invariance as is discussed for k~ in the previous section. Note that we obtain inter- 
action terms through this procedure. However as we will see just below, the forms 
of their continuum limit are completely determined by the dimension, parity and 
5*0(8) invariance. First we concentrate on free parts of the supercharges Q^^^ and 
Q^'^\ i.e., ignore the interaction terms. By using the power counting, Eqs. (2-39) and 
(2-47), SO(8) invariance and the parity symmetry, Eqs. (2-42), (2-43) and (2-48), 
we can deduce the following forms of free supercharges in the e ^ limit: 

QilL = J daa,ea{a), 

QirLa = j daa2Xa{cr), 

Q^L = J da{b2k\a)jlXa{a) + C2x'\aHjaia)). (2-50) 



18 H. Aoki, S. Iso, H. Kawai, Y. Kitazawa, A. Tsuchiya and T. Tada 

In Eq. (2-50) we have excluded terms such as ^x*7*A by translation invariance. It 
is easy to see that all possible terms which appear with negative powers of e are 
forbidden by the symmetries. In this sense the existence of the continuum limit 
is guaranteed by the symmetries. We can also fix undetermined coefficients in Eq. 
(2-50) by the M=2 supersymmetry (2-49) as follows. From {Qa\Q^^^} = V^P'^dab, 
} = P^T^d and qI'^} = P^a, we obtain 

aia2 = V2, aib2 = 1 and a2Ci = 1. (2'51) 
Therefore Eq. (2-50) is reduced to 

<5freea = «l / d<70ai<7), 

Q^a = ^ / daXaia), 

QSL = jdcy [h,\aH,\a{a) + C2x'\a)^Mcy)) • (2-52) 

The free part of the Hamiltonian / dak~{a) is obtained by {Q^(^\Q^^^} = V^H5^^ 
as 

Iffree = jda {^k\af + l=h,C2x'\af - i^\'^{a)\a{a) - ^^c2eM0a{a)) . 

(2-53) 

In order to compare these results with the Green-Schwarz light-cone formalism, we 
redefine the fermionic variables as 

>^a = x[^r]Xa and e„ = W-7?*^a, (2-54) 

where 77 = e^. We also introduce rescaled supercharges and by 

Q(i) = y^r/*Q^ and Q^^^ = -^Q^^\ (2-55) 

V«iOi 

In terms of these new quantities, Eqs. (2-52) and (2-53) become 

Qfreea = j da0a{cr), 

CL = daUa), 

Qieea = j d^ (^-^k\a)-liJa{(T) + ix\a)-ilX{(^)) , 

dd = / da{ieiaUX{^)-^hc2x'\aUX{a)), 

(2-56) 
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which completely agree with the light-cone Green-Schwarz free Hamiltonian and su- 
percharges for type IIB superstring. This fact also justifies the analytic continuation 
introduced for fcrmionic fields in Ref. 1). Wc also note that we have obtained the re- 
lation 61 C2 and 61C2 should be equal to e/g^I multiplied by some numerical 
constant as is illustrated in the previous subsection. 

2.4.4. Interaction parts of supercharges and Hamiltonian 

In this subsection, we examine the structure of the interaction parts of the super- 
charges and the Hamiltonian. First we consider the contributions of order 1//, which 
correspond to 3-Reggeon vertices in string field theory. Since our free Hamiltonian 
is equal to that of the Green-Schwarz light-cone formalism and the interactions of 
loops are local in our loop equations, we can use the same arguments as in light-cone 
string field theory. In general, the operators inserted near the interaction points 
in 3-Reggeon vertices generate divergences coming from the Mandelstam mapping. 
Since our Wilson loops are written by the variables fc* and A, the corresponding 3- 
Reggeon vertices should consist of delta functions representing the matching of three 
strings in the k — \ space, which is the same as in Ref. 18). Therefore the fc*, x'* 
and Xa diverge as 1 /^/e near the interaction points while 9a is of order e'^ there. We 
also note that every derivative of a acting on the fields introduces an extra factor of 
l/s/e. Therefore the interaction part at order 1/7 of the supercharges possesses the 
following general structure: 

J J da J da,e\kT{x'Y{XanOaY 

X (derivative) '^(products of delta functions for A;' and Xa), (2-57) 

where A;*, x'*, Xa and 6a represent the operators inserted near the interaction points, 
and K is the total number of derivatives acting on these operators. 

For example, let us consider the interaction part of Q^^^ ■ In this case, the 

dimensional analysis [Qint^] = [Qfreedl leads to -rj + a + P + ^'j + ^6 + k - 1 = ^, 
where we used the relation [/] = —1 and the fact that the dimensions of the delta 
functions are canceled due to supersymmetry. Therefore the total powers of e which 
appear in the interaction part of Q^^^ is evaluated as 

The case in which a = f3 = 5 = K = is excluded by S0{8) invariance. We can 
consider four cases in which ( = —1: {l)a = 1 and P = 5 = k = 0, {2)13 = 1 and 
a = (5 = K = 0, (3)(5 = 1 and a = /3 = k = and (4)k = 1 and a = (3 = 5 = Cases 
(3) and (4) are not permitted by 5*0 (8) invariance. If we take the large- iV limit with 
le kept fixed. Cases (1) and (2) survive in the e — > limit. This limit corresponds to 
taking g^t ~ l//e. Note that in this limit all of the other cases vanish because ^ is 
larger than —1 for them. Furthermore we can restrict the values of 7 by the parity 
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symmetry and deduce the structure of Qintd ^ fohows: 

Qi2a = tJ'^''J da,{^/-ek\{V~eXa? + (VeK)') + Vex'^VeXa + {V~eKf)) 
X (products of delta functions for fc* and Xa)- (2-59) 
This structure agrees with that of the light-cone string field theory. Applying a 

(2) 

similar analysis to Qint a; we obtain 

QSd = tJ'^''J da{V~ek\V~eXa + {^feXaf) + Vex"{{V~eXaf + (VeXaf)) 
X (products of delta functions for fc* and Xa), (2-60) 

which also agrees with the light-cone string field theory. As for Q^^ and Q^a \ no 
l/I contribution remains non-zero in this limit since the minimum value of is — ^ 

in these cases. Therefore we conclude that Qj^ta ^^'^ Qinta equal to zero at 
order l/I, which is again consistent with the light-cone string field theory. Note 
that the right-hand sides of Eqs. (2-59) and (2-60) are uniquely determined by M=2 
supersymmetry, as is shown in Ref. 18). Finally the anti-commutation relation 
{Q^^\Q^^^} = ^HS^j^ fixes the interaction part of H, which is certainly consistent 
with the light-cone string field theory. 

Next we consider the contributions of order l/I^ [k > 2), which correspond to 
(A; -I- 2)-Reggeon vertices. The general structure of the interaction part is represented 
as 

j^jdajdai- ■■dake'^{kT{x'Y{Xa)nea)'{deTiYatiYe)'' 

X (products of delta functions for fe* and Xa). (2-61) 

From the Mandelstam mapping in these cases, it is natural to consider that the k^, 
fc_ 

X and Xa diverge as e fe+i near the interaction points. Therefore the total power ^ 
of e is evaluated as 

where p is ^ for Qi2,a Qinta' 1 ^^'^ Qinta Qinta ^^'^ terms in which 
C < —k survive in the e — limit if le is fixed. It is verified easily that there are 
no surviving terms for any values of k in Qinta and Q^f^-a e — > limit, which 

is consistent with the light-cone string field theory. Using 50(8) invariance, we can 
show that in Q^^^ ^ and Q^^^ ^ some terms with 7 equal to five might survive for = 2 
and ones with 7 equal to seven for A; = 2 and k = 3. Presumably it is not possible 
to satisfy M=2 supersymmetry only by these restricted terms. Therefore we may 
conclude that there are no contributions of order l/I^ {k > 2) in Qintd' Qinta 
the Hamiltonian, which is also consistent with the light-cone string field theory. 

In this way, we almost confirm that our IIB matrix model reproduces the light- 
cone string field theory for type IIB superstring. 
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We have the following two relations about the scaling limit. 

9^1 



'2 

a ~ 



(2-63) 



(2-64) 



Since we assume that / ~ A'^" and e ^ N ^, we obtain from Eq. (2-63) the prescrip- 
tion of the scaling limit: 

g'^N"-^^ ~ const. (2-65) 
Prom Eq. (2-64), in order to obtain a finite string coupling, we also have a relation 

a = b. (2-66) 

We need to estimate large N behavior of / and e, which should be consistent with 
Eq. (2-66), in order to fix the prescription of the scaling limit completely. Assuming, 
for example, that e ^ \fa' /R and R ~ ^N^, we obtain as a candidate 



9 i 

g Ni ^ const. 



§3. Dynamics of eigenvalues and space-time generation 



(2-67) 



In this section we analyze the structure of space-time, and in particular, try 
to explain why our space-time is four-dimensional. As we mentioned in the in- 
troduction, the diagonal elements of the bosonic matrices can be interpreted as 
space-time itself. For example, if the diagonal elements distribute within a manifold 
which extends in four dimensions but shrinks in six dimensions, then a natural in- 
terpretation is that the space-time is four-dimensional. We thus derive an effective 
theory for the diagonal elements, and analyze their distribution. 

We decompose into diagonal part and off-diagonal part A^. We also 
decompose ?/> into diagonal part ^ and off-diagonal part il): 



( e 



I X 



\ 



X, 



\ 



< ) 



V 



(3-1) 



where and satisfy the constraints a;^ = and = 0, respectively, 

since we may fix the U{1) part by translation invariance. We then integrate out the 
off-diagonal parts A^ and '0 and obtain the effective action for supercoordinates of 
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space-time SefffX, ^]. The effective action for the space-time coordinates 5'eff[X] can 
be obtained by further integrating out ^ : 

= J dXe-^<=ff[^l, (3-2) 

where dX and stand for Jlil:!^ 11^=0 '^^fj T[f=i^ Il^=i ^Cai respectively. 

We perform integrations over off-diagonal parts A^ and ip by the perturbative 
expansion in g"^, which is valid when all of the diagonal elements are widely separated 
from one another: |x' — S> g^^^- After adding a gauge fixing and the Faddeev- 
Popov ghost term 

-5g.f. + Sf.p. = -^Tri[X^,A^]^) - ^Tri[X^,b][A^,c]), (3-3) 

the action can be expanded up to the second order of the off-diagonal elements A^, ijj 
as 

HC - i')r^i>''Af + Tr^iC - • (3-4) 

The first and the second terms are the kinetic terms for A and V' respectively, while 
the last two terms are A-ij)^ vertices. A bosonic off-diagonal element is transmuted 
to a fermionic off-diagonal element ■0*"' emitting a fermion zeromode ^* or This 
vertex conserves indices for space-time points, i, j. Note that the propagators for 
An and ■0 behave as r7^i,/(x* — x^Y (^Jx ~ ~ respectively, thus 

they decrease in the long distances. 

We obtain the effective action for the zeromodes, and at one- loop level: 



/ 



dAd-^^dhdc e-(^2+5g.f.+5F.p.) = JJ det^^ (ry^'' + S'^^.^ 

i<j 

^e-^ey°°^[^'«l, (3-5) 



where 

5g^)=^^T--r^. (3-6) 

Here ^^-^ and are abbreviations for ^* — and — x^. The effective action can 
be expanded as 

i<j 

= + (3-7) 
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which is a sum of all pairs {ij) of space-time points. Here the symbol tr in the lower 
case stands for the trace over Lorentz indices, ii, v. Other terms in the expansion 
vanish due to the properties of Majorana-Weyl fcrmions in ten dimensions. 

Note that the one- loop effective action S^q°°'^[X,^] has J\f = 2 supersymmetry, 

f^Wx^ =ie-ir^e f5(2)x^ =0 

U«e =0 ' =62 ' ^ ' 

which is a remnant of the one in the original theory, 

jdWA^ =ieir^i; r5(2)A^ =0 q^ 

Transformations for ei = £2 and ei = —€2 correspond to those generated by A/" = 1 
supersymmetry generator Q and its covariant derivative D. In this sense zeromodes 
of Xi and may be viewed as supercoordinates of = 1 superspace. 

The effective action for the space-time coordinates is given by further inte- 
grating out the fermion zeromodes ^f. 



= I dXd^H 



(3-10) 



Here the products are taken over all possible different pairs (ij) of the space-time 
points. When we expand the multi-products, we select one of the three different 
factors, 1, tr(5|'.^.))/4 or {tT{Sf..^)/8 + (tr(5j^.)))V32) for each pair of (ij). Since 

the last two factors are functions of (x^ — xjj), they can be visualized by bonds 
that connect the space-time points x^ and x^. Since the factors tr{S^-j-^)/A and 
(tr(S'^^^.-|)/8 -|- (tr(S'^.^p)^/32) contain 8 and 16 spinor components of , we call 
them an 8-fold bond and a 16-fold bond, respectively. We do not assign any bond 
to the factor 1. In this way we can associate each term in the expansion of multi- 
products in Eq. (3-10) with a graph connecting the space-time points by 8- fold bonds 
and 16-fold bonds. Therefore the multi-products in Eq. (3-10) can be replaced by a 
summation over all possible graphs: 



/ 



dxe-^l^'°°'m = J dxdc E n 



G;graph (ij):bond of G 

' '<^y ] /i (^^jsyv , tr(^s,))' 



(3-11) 

Here we sum over all possible graphs consisting of 8-fold and 16-fold bonds. For each 
bond (ij) of G, we assign the first or the second factor depending on whether it is 
an 8-fold or a 16-fold bond. 
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In order to saturate the grassmann integration d^, we need 16(A'' — 1) fermion 
zeromodes. Since an 8-fold bond and a 16-fold bond contain 8 and 16 fermion 
zeromodes respectively, graphs remaining after the integrations are those where 
the sum of the number of 8-fold bonds and twice the number of 16-fold bonds is 
equal to 2(A'^ — 1). Thus, the number of bonds is of order N, which is much smaller 
than possible number of pairs N{N — l)/2. In this sense N space-time points are 
weakly bound. 

Since 8-fold bond and 16-fold bond terms behave as {x^^)~^'^ and (a;*-')^^'^, both 
are strong attractive interactions. Hence, only the closer points can be connected 
by the bonds, and these interactions become local. On the other hand, since all 
possible graphs must be summed up, this system has a permutation invariance among 
the points. This reminds us of summation over all triangulations in the dynamical 
triangulation approach to quantum gravity. We come back to this analogy in §4. 

In order to see some important features of the system (3-11), let us take an 
approximation. If there were only 16-fold bonds, considerable simplifications take 
place and the ^ integrations can be performed exactly. Since a 16-fold bond term 
contains 16 fermion zeromodes, it is proportional to delta function of the grassmann 
variables as 

tr(4,))~5^^«)(r-a^-^-pj. (3-12) 

If there is a loop in the graph, the contribution vanishes, since the product of delta 
functions of grassmann variables on the loop vanishes: 

Also, as we mentioned above, the remaining graphs have {N — 1) number of 16-fold 
bonds. Hence, the remaining graphs arc tree graphs which connect all N points. 
Such type of graphs are called "maximal trees" . We also note that all maximal trees 
contribute equally as we can see by performing ^ integrations from the end points 
of each maximal tree. Therefore, only with 16-fold bonds, the distribution of 
becomes "branched polymer" type: 

= E n SHiWj- (3-14) 

Grmaximal graph (ij):bond of G ^ ^ 

Note that all points are connected by the bonds, and each x^^ integration can 
be performed independently and converges for large x^^ on each bond. Thus this 
system is infrared convergent. As we proved rigorously in Ref. 3), this feature of IR 
convergence holds even with 8-fold bonds, and also, to all orders in perturbation ex- 
pansion. This is consistent with the explicit calculations of the partition function. 
This shows that all points are gathered as a single bunch and hence space-time is in- 
separable. Thus, the size and the dimensionality of the space-time can be determined 
dynamically. Note also that the dynamics of branched polymer is well known and 
its Hausdorff' dimension is four. It is conceivable that the smooth four-dimensional 
space-time emerges by taking account of the effects from 8-fold bonds as we argue 
in what follows. Therefore, the model (3- 11) constitutes a candidate of models for 
dynamical generation of four-dimensional space-time. 



/ 
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Before going into the analysis of the space-time structures by using the effective 
action (3-11), which is vahd in the long distance, we consider the short distance 
behavior of the system. Let us suppose that a pair of the bosonic coordinates are 
degenerate but the rest of the coordinates are well separated from one another and 
from the center of mass coordinates of the pair. We can determine the dynamics 
of the relative coordinates of the pair of the points, from the exact solution for the 
SU (2) case. The distribution for the relative coordinates is 




<^g ' ^ ' 

We conclude that there is a pairwise repulsive potential of —8 Inr type when two 
coordinates are close to each other. It is clear that these considerations are valid 
for arbitrary numbers of degenerate pairs although the center of mass coordinates 
should be well separated. Although it is possible to repeat these considerations to 
the cases with higher degeneracy, the analysis becomes more complicated. Therefore 
we choose to adopt a phenomenological approach and assume the existence of the 
hardcore repulsive potential of the following form: 

5core[^] =E5(^'-^'), (3-16) 

where 

g{x -x)-<^^ for {x' - x^ f » g ' ^^'^^^ 

Hereafter, we investigate the structure of space-time by using the one loop effective 
action (3-11) plus the phenomenological hardcore potential (3-16). 

If the number of 16-fold bonds is much larger than that of 8-fold bonds, the 
system behaves as a branched polymer. A small number of 8-fold bonds may fold the 
branched polymer into a lower-dimensional manifold. This can happen as in protein, 
a chain of amino acids is folded into a lower-dimensional object like /3-sheet, by 
perturbative interactions. Since the Hausdorff dimension of the branched polymer is 
four, the core interactions exclude the manifolds in less than four dimensions. Thus, 
four-dimensional space-time can be realized by this mechanism. We are checking 
this conjecture by numerical simulations, which we mention later in this section. 

On the other hand, if the 8-fold bonds dominate, the number of bonds are of 
order 2A'^, twice as much as in the above case. Thus, the entropy of graph rear- 
rangement becomes more important, and the system might behave as a mean field 
phase, where all the N points condensate into a finite volume. However, the core 
interactions prohibit an infinite density state, and the system behaves as a droplet. 
The 8-fold bond interaction can be written as 

tr(%))4 cx C>^^^P-^-'^s^l . . . elV%J{.x^^f\ (3-18) 
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where C'^^^^^i-^s is an invariant tensor, and V^'l^^ is a fourth rank symmetric 
traceless tensor constructed from x^. Since V^^;^p is traceless, the average over 
orientations of gives a suppression factor for each 8-fold bond. This suppression 
factor becomes weak if the system becomes lower-dimensional. Also, the inte- 
grations give contractions of V^j^;^p among different bonds, which consist of inner 
products x^^x^^ between different bonds. These angle-dependent interactions may 
favor lower-dimensional space-time. 

We can study self-consistently which phase is realized; 16-fold bond dominant 
branched polymer phase or 8-fold bond dominant droplet phase. The phase with the 
lowest free energy is realized. However, in any phase, four-dimensional space-time 
can be realized by one of the above mentioned mechanisms or by some combinations 
of them. 

In the remainder of this section, we show how we perform numerical simulations. 
Our conjecture is that a small number of 8-fold bonds fold the branched polymer 



2.5 




' • ' • ' • ' • ' 

50 100 150 200 

iterations/1000 

Fig. 2. Result of the numerical simulation, where = 800 and the number of 8-fold bond is fixed 
to be 140. The horizontal axis represents the number of iterations, the vertical axis the length 
of the ten principal axes in ten dimensions. Two of them are relatively large, suggesting the 
existence of anisotropy in space-time. 
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into four-dimensional space-time. Thus we take the following model: 

„ Af-l 

^ = / n d''^' E E (3-19) 

j=l G:maximal tree B:sets of 8-fold bonds 

where 

S= 121n[(a:^-a;^)2 + l]+ ^ Q\u[{x^ - x^f + I] 

+ Y.e{l- (x^^y)(-4) H{x^^f{2 - {x^^f))]. (3-20) 

Here we fix the number of 8-fold bonds by hand. It is enough to check the above 
mentioned conjecture, although the number of 8-fold bonds is actually fixed by dy- 
namics. 

Wc generate distribution of {a;^} by the Monte Carlo method. We then measure 
moment of inertia, 

N 
i=l 

and diagonalize it. The ten eigenvalues are the length squared in the principal axes 
in ten dimensions. In this way, we analyze anisotropy of tcn-dimcnsional space-time. 
For example, if four of the eigenvalues are much larger than the others, and grow as 
we take TV large, it means space-time is four-dimensional. 

Figure 2 is our preliminary result. In AT = 800, two of the ten eigenvalues 
are relatively large, suggesting the existence of anisotropy. We hope to see four 
eigenvalues become larger as we increase N. However, in order to see this, we need 
at least 4^-2^ ~ 16, 000 points, which is quite a large number in the current computer 
power. Thus, we arc trying to make some modifications to the model (3-20) to study 
the system with larger N. The results will be reported soon. 

§4. Symmetries in the low energy theory 

4.1. Local gauge invariance 

Once we describe the space-time as a dynamically generated distribution of the 
eigenvalues, low-energy effective theory in the space-time can be obtained by solving 
the dynamics of the fluctuations around the background X^. Both of the space-time 
Xfj^ and matter are unified in the same matrices and should be determined 
dynamically. Low-energy fluctuations are generally composites of and V') and it 
is natural from the analysis of the Schwinger-Dyson equation for the Wilson loops 
that a local operator in space-time is given by a microscopic limit of the Wilson loop 
operators, such as 

w{k;0) = Tr[0(^,V)exp(iA;^A^)]. (4-1) 

Here, 0{A,ip) is some operator made of A^ and ip. In order to identify the total 
momentum of this operator with k, the operator 0{A, ijj) should be invariant under 
a constant shift of A^, that is, a translation in the space-time coordinates. 
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In the first approximation around the diagonal background X^, the coordinate 
representation of this operator is given by 



w 



{x;0) = I 



(27r) 



10 



exp(— zfe'^x^) w{k; O) 



N 



(4-2) 



i=l 



Here, we have replaced A^ by + ^4^ and taken the leading term. On is the {ii) 
component of the operator O. Due to the delta function, the operator has support 
only in the domain where the eigenvalues distribute. Vanishing of the operator 
w{x; O) outside of the domain of the distributed eigenvalues implies that space-time 
simply does not exist outside the domain. This fact supports our interpretation of 
the space-time in IIB matrix model. 

We can apply a similar analysis to strings which propagate in the space-time. In 
the expansion, the correlation between Wilson loop operators can be evaluated 
by summing over all surfaces made of Feynman diagrams connecting the Wilson 
loops at the boundary. This surface is interpreted as a string world sheet connecting 
strings at the boundaries. Each eigenvalue (xj^) associated to a loop in the diagrams 
represents a coordinate on the world sheet, and it takes a value in the eigenvalue dis- 
tribution in the leading approximation around the diagonal background X^. Hence, 
a string world sheet evolves only in the space-time of the eigenvalue distribution, 
which again supports our interpretation of space-time. 

It is generally difficult to obtain how fluctuations propagate in the eigenvalue 
distribution, which is reminiscent of the QCD effective theory: In QCD excitations 
are expressed as composite operators of microscopic variables, and their low-energy 
dynamics can be discussed only through a symmetry argument, namely an argument 
based on chiral symmetry. Also, in our case, we will show that there are eigenvalue 
distributions around which symmetry arguments allow us to discuss the low-energy 
dynamics for some excitations. Suppose that the eigenvalue distribution forms clus- 
ters consisting of n eigenvalues. At a length scale much larger than the size of each 
cluster, the SU{N) symmetry is broken down to SU{n)'^, where m = N/n. 

We can expand A^ and i{j around such a background similarly to the analysis 
in the previous section. First, write A^ and in block forms: 

/ Af Af ... \ 

A^} 



A,,= 



V 

^21 ^22 



V 



Aram 



(4-3) 



Each block, A^^ or ■0'-' , is an n x n matrix, and the diagonal blocks can be further 
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decomposed: 

^n=^il+^i\ (4-4) 

where 1 is an n x n unit matrix and tr = 0. Here, tr means the trace for a 
submatrix of n x n. We interpret each chister of the eigenvalues as being a space- 
time point with an internal structure SU{n). Since each SU{n) symmetry acts on 
the variables at position i independently, the unbroken SU{n)"^ symmetry can be 
regarded as being local gauge symmetry. Indeed, under a gauge transformation of 
the unbroken SU{n)"^ symmetry, 

/ gi \ 



92 



V 9m J 



G SU{nr C SU{N), (4-5) 



the diagonal block fields, A^^ and transform as adjoint matters, (i.e., site variables 
in the lattice gauge theory), while the off-diagonal block fields, A'j^ and transform 
as gauge connections (i.e., link variables): 

r^9ir9}, 

A^^9iA^9l, 

4^'' ^ 9i4^''9]- (4-6) 

Some of the dynamics for low-energy excitations is governed by this local gauge 
invariance. Gauge fields live on the links and transform as the link variables in lattice 
gauge theory. In our case, we have too many such fields (at least 10 boson fields A^^ 
for a link (ij) ), but only one unitary link variable is assured to be massless by the 
gauge symmetry, and the others acquire mass dynamically. Therefore, in deriving 
low-energy effective theory, we first apply a polar decomposition to into unitary 
and hermitian degrees of freedom, and identify all of the unitary components of 
by setting them to be one common field on each link. Wc have to integrate those 
massive off-diagonal block fields A^^ while keeping the unitary components and 
this procedure is performed by the following replacements: 

^ = 0, 

A^^A^ = 9'j^U'' ^ U^'- (4-7) 

In the second equation, the factor, l/(x*-')^, corresponds to the propagation of the 
hermitian degrees of freedom while the appearance of the link variable C/*^ corre- 
sponds to keeping the unitary degrees of freedom. Higher order correlations can be 
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obtained by using Wick theorem in general except ten-body correlation function for 
A]^ . Due to the chiralness of the ten-dimensional fermion, we obtain an extra term 
proportional to e^^...^!,,, coming from fermion one-loop integral. 

In order to derive the effective theory for fluctuations around the assumed back- 
ground, we integrate massive fields first and obtain effective action for other fields 
as we have obtained the effective action for the diagonal components (i.e., space- 
time coordinates) in the previous section. Generally speaking, we can expect any 
terms which are not forbidden by symmetries; supersymmetry and local SU (n) gauge 
symmetry. 

A plaquette action for gauge fields U^^ can be generated as follows. The relevant 
terms in the action are 

54 = ^ ^ tT{Ai:iAiJ^-A^Al!:)A^A^. (4-8) 

By integrating out the hermitian degrees of freedom of the off-diagonal blocks with 
the procedure (4-7), this action vanishes; S'4 = 0. However, interactions generated 
by {84)^ induce a kinetic term for the gauge field: 

- ^^E^^ J^^IJ^T^J^f^^ (4-9) 

This is the plaquette action generated by a Wilson loop for the adjoint representation, 
and hence the gauge field U^^ indeed propagates in the space-time of the eigenvalue 
distribution. The gauge field can hop between any pair of space-time points, but the 
hopping is suppressed by for distant points and we will recover locality in the 
continuum limit. Similarly we can obtain a gauge invariant hopping term for adjoint 
fermion 

To summarize this subsection, supposing that the distribution of the eigenvalues 
consists of small clusters of size n, we have shown that the low-energy effective theory 
contains several massless fields, such as the gauge field associated with the local 
SU{n) gauge symmetry and fermion field in the adjoint representation of SU (n) 
gauge symmetry. Gauge-invariant kinetic terms were also derived. Our findings 
are reminiscent of those which are obtained by considering n coincident D9-branes. 
Presumably our argument here is related to the standard argument of coincident 
D-branes. 

Our system is a lattice gauge theory on a dynamically generated random lattice. 
It is invariant under a permutation for the set of the m discrete space-time points, 
since the permutation group Sm is a subgroup of the original SU (N) symmetry. The 
existence of the permutation symmetry is the crucial difference from the ordinary 
lattice gauge theory on a fixed lattice, which becomes important in deriving the 
diffeomorphism invar iance of our model. We will come back to this point in the next 
subsection. Although the permutation invariance requires that all space-time points 
are equivalent, locality in the space-time will be assured due to suppression of the 
hopping term between distant points. In general, however, we need a sufficient power 
for the damping of the hopping terms in order to assure locality in the continuum 
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limit. Though we do not yet know the real condition for locality, we expect that 
terms with lower powers are canceled due to supersymmetry or by averaging over 
gauge fields. 

4.2. Diffeomorphism invariance 

As shown in §2, the one-loop effective action for the space-time points is de- 
scribed as a statistical system of N points whose coordinates are x^. Integration 
over the fermion zeromodes ^ gives the Boltzmann weight, which depends on a graph 
(or network) connecting the space-time points locally by the bond interactions: 

Z= J2 JdXW[X;G]. (4-10) 

G;graph 

VF[X;G] is a complicated function of a configuration X and a graph G. An im- 
portant property is that the weight is suppressed at least by a damping factor of 
l/(x* — x^)^^, when two points, i and j, are connected. This system is, of course, 
invariant under permutations S^*^ of N space-time points, which is a subgroup of 
the original symmetry SU{N), while the Boltzmann weight for each graph G is not. 
The invariance is realized by summing over all possible graphs. In other words, 
the system becomes permutation invariant by rearrangements of the bonds in the 
network of the space-time points. This reminds us of the dynamical triangulation 
approach to quantum gravity, "^^^ where diffeomorphism invariance is believed to arise 
from summing all possible triangulations. Our system satisfies both the locality and 
permutation invariance simultaneously by summing over all possible graphs whose 
points arc connected through the local interactions. 

Now we see that the permutation invariance of our system actually leads to 
diffeomorphism invariance. To see how the background metric is encoded in the 
effective action for low-energy excitations, let us consider, as an example, a scalar 
field 0* propagating in distributed eigenvalues. The effective action will be given by 

5 = E ^^^^/(^^ - + E "^('^^)^ (4-11) 

where f{x) is a function decreasing sufficiently fast at infinity to assure locality in 
the space-time. Introducing the density function of the eigenvalues, 

(x — X*), (4-12) 

i 

and a field (/)(x) which satisfies (/)(x*) = the action can be rewritten as 

S = J dxdy{p{x)p{y)) ('^(^)-'^(^))' /(^ -y) + mj dx{p{x))<l>{xf . (4-13) 

In this section we consider general eigenvalue distributions in which all eigenvalues have 

nondegenerate space-time coordinates. If we take the cluster type distribution considered in the 
previous subsection, the permutation symmetry responsible for the diffeomorphism invariance should 

be ySjv/n- 
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Here, the expectation, (•••), for the density and the density correlation means that we 
have taken average over configurations X and networks G of the space-time points. 
Normalizing the density correlation in terms of the density, 

{p{x)p{y)) = {p{x)){p{y))r{x,y), (4-14) 

and expanding (f>{x) — (j){y) = {x — y)^d'-^4'{x) + • • •, the action becomes 



^ = \j dx{p{x)) J dy{p{y)){x - y)^{x - y)vf{x - y)r{x, y) 



d^(t){x)d''(t>{x) 



+m 



J dx{p{x))(t){xf (4-15) 



This expansion shows that the field 4){x) propagating in the eigenvalue distribution 
feels the density correlation as the background metric, while the density itself as 
vacuum expectation value of the dilaton field. Namely, we can identify 



9,iu{x) J dy{p{y)){x - y)f,{x - y)i,f{x - y)r{x,y), (4-16) 

Vde-''^'^ - {p{x)). (4-17) 

If the density correlation respects the original translational and rotational symmetry, 
that is, if they are not spontaneously broken, the metric becomes flat, Qfj^u ~ rj^i,. 
(Normalization can be absorbed by the dilaton vev.) The fact that the background 
metric is encoded in the density correlations indicates that our system is general co- 
variant, even though the JIB matrix model action (1-1) defined in flat ten dimensions 
does not have a manifest general covariance. 

Then, let us see how the diffeomorphism invariance is realized in our model. 
The action (1-1) is invariant under the permutation Sn of the eigenvalues, which is 
a subgroup of SU{N). Under a permutation, 

x^^x^W for aeSN, (4-18) 

the field transforms into (/)'^(*). Then, from the definition of the field (f>{x), we 
should extend the transformation (4T8) into x, 

X ^ (4-19) 

such that ^(x*) = x"^^\ Under this transformation, the eigenvalue density trans- 
forms scalar density and the field (p{ ctS cl scalar field. On the other hand, 
the metric transforms as a second-rank tensor, if the function /(x — y) decreases 
rapidly around x = y and the y integral in Eq. (4-16) has support only near y = x. 
The tensor property of the metric is also required from the invariance of the action 
under transformation (4-19). In this way, the invariance under a permutation of the 
eigenvalues leads to the invariance of the low-energy effective action under general 
coordinate transformations. 

The background metric is encoded in the density correlation of the eigenvalues. 
Since we have started from the IIB matrix model action (1-1) which is Poincare- 
invariant, the density correlation is expected to be translational and rotational in- 
variant, and we may obtain a low-energy effective action in a fiat background. A 
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nontrivial background can be induced dynamically if the Lorentz symmetry is spon- 
taneously broken and the eigenvalues are nontrivially distributed. 

A nontrivial background can also be described by condensing a graviton opera- 
tor. Bosonic parts of graviton and dilaton operators are given by 

S^,{k) ~ IV(F^;,F\e^^-^) + v\ (4-20) 
D{k) ~ Tr(F2e^'=-^). (4-21) 

Their condensation induces extra terms in the IIB matrix model action, 

-Scond = Jdk(j2 h>'''{k)S^Ak) + h{k)D{k)) . (4-22) 

We can similarly obtain an effective action for fluctuations around a diagonal back- 
ground from this modified matrix model action. Condensation of dilaton changes the 
Yang-Mills coupling constant g locally in space-time. Since g is the only dimension- 
ful constant in our model, and thus determines the fundamental length scale, a local 
change in g leads to a local change in the eigenvalue density. This is consistent with 
our earlier discussion that the dilaton expectation value is encoded in the eigenvalue 
density. On the other hand, the condensation of graviton induces an asymmetry of 
space-time. For a condensation of the A; = graviton mode, it is obvious that the 
condensation can be compensated by a field redefinition of matrices 

A^^{6; + h''^)A,, (4-23) 

and the two models, the original IIB matrix model and the modified one with the 
k = graviton condensation, are directly related through the above field redefini- 
tion. The density of the eigenvalues is mapped accordingly, and the density cor- 
relation is expected to become asymmetric in the modified matrix model. For a 
more general condensation, if the graviton operator S^y{x) (coordinate representa- 
tion of Eq. (4-20)) changes only the local property of the dynamics of the eigenvalues, 
the density correlation will become asymmetric locally in space-time around x, and 
therefore induces a local change in the background metric. 

Our low-energy effective action is formulated as a lattice gauge theory on a dy- 
namically generated random lattice. Since the lattice itself is generated dynamically 
from matrices, we must sum over all possible graphs. In this way, our system is 
permutation Sm invariant, which is responsible for the diffeomorphism invariance. 
The background metric is encoded in the density correlation of the eigenvalues, and 
the low-energy effective action becomes manifestly general covariant. The graviton 
operator is represented as fluctuation around the background space-time, and is con- 
structed from the off-diagonal components of the matrices. A microscopic derivation 
of the propagation of the graviton is difficult to obtain, but once we have clarified 
the underlying diffeomorphism symmetry, it is natural that the low-energy effective 
action for the graviton is described by the Einstein Hilbert action. By employing 
this diffeomorphism invariance and the supcrsymmetry, we will be able to derive the 
low-energy behavior of the graviton multiplets, which will be reported in a separate 
paper. 
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§5. Discussion 



We have reviewed the current status of the type IIB matrix model, which is 
proposed as a constructive definition of supcrstring. 

There are still several conceptual issues. We have obtained a nonabelian gauge 
symmetry from the IIB matrix model by assuming a particular eigenvalue distribu- 
tion. This indicates that this vacuum is not a perturbative vacuum of the type IIB 
superstring. Instead, we may wonder if this is a perturbative vacuum of a hctcrotic 
string or a type I string realized in a nonperturbative way within the IIB matrix 
model. As we discussed in the introduction, our matrix model contains both of 
the world sheets of the fundamental IIB string and those of the D-strings. By a 
semiclassical correspondence (1-8), we have identified a IIB supcrstring in the Schild 
gauge where tr is interpreted as integration over a D-string world sheet. We can also 
construct an F-string world sheet in terms of surfaces made of Feynman diagrams, 
whose SU{N) index represents the space-time coordinate of a world sheet point. In 
both cases, if we assume an eigenvalue distribution consisting of small clusters, an 
internal structure appears on the world sheet, and hence current algebra may arise 
and there is a possibility to describe heterotic string within the type IIB matrix 
model. 

Another issue is how to describe global topology in the IIB matrix model. A 
simple example is a torus compactification. A possible procedure of torus compact- 
ification^^) is to identify with A^ + by embedding a derivative operator into 
our matrix configuration. Therefore, is taken as infinity from the beginning. Since 
this procedure has a subtlety in the large N limit, we require a careful examination 
of the double scaling limit. 

We also do not yet know how we can describe chiral fermions in lower dimensions 
after compactification. If we naively consider a low energy effective theory on a 
four-dimensional space-time generated by distributed eigenvalues, we will obtain 
fermions with both chiralities. A possible mechanism to produce a chiral fermion is 
to consider a compactified six-dimensional space-time with a nontrivial index, or a 
parity violating background. We should replace the n x n block matrices considered 
in §4 by (n -|- oo) x {n + oo) matrices. Size n part of a block represents the low energy 
gauge symmetry as before. The rest of the infinite size represents an internal space 
with six dimensions, which should have a non-trivial index. For simplicity, let us 
consider a two-dimensional internal space. Then the following background 



gives the {ii) block components of a desired background. Here a = 1 ~ 4 and P, Q 
are infinite dimensional matrices satisfying [P, Q] = —i. This background is invariant 
under U{n) x U{1) in each block and this becomes the local gauge symmetry in four- 
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dimensional space-time. The effective theory for low energy fluctuations around the 
background are similarly obtained as in §4. All the off-diagonal fields become massive 
and we integrate them over except a gauge field. In diagonal blocks, there arc several 
fields that can be massless. Writing fluctuations in a diagonal block, O = AJ^ or 

as 



o = rx^ , (5-1) 




Oa transforms as an adjoint representation, Oj^ as a fundamental representation and 
Os is a singlet for SU{n) transformation. P or Q acts like PO^p on O]^, and [P, 0$] 
on S-components. Then we can show that we obtain a massless chiral fermion 
with fundamental representation for gauge symmetry SU{n), whose wavefunction 
in (P, Q) space is given by the groundstate wavefunction of a harmonic oscillator. 
Other possibly massless fcrmions arc vector-like and will acquire mass unless they are 
protected by supersymmetry. This is the simplest way to obtain a chiral fermion in 
four-dimensional space-time. In this construction, all non-singlet fields transforming 
as adjoint or fundamental representations live in four dimensions and localized at 
= = 0. However, singlet fields including graviton propagate in the bulk (here, 
six dimensions). In order to have a four-dimensional theory, we need to compactify 
the internal six-dimensional space. This is discussed in a separate paper. 

It is also desirable to construct AdS type backgrounds in our approach. Let 
us recall the metric of AdS^ x S^: 

ds'^ = ^{dx^ + dz^) + R^dQl. (5-2) 

The volume factor yfg ~ 4^ is sharply peaked at z = in AdS^. Since we have 
argued that ^fg is proportional to the density distribution in IIB matrix model, 
such a background may be represented by the eigenvalue distribution which is also 
sharply peaked at 2; = namely at the four-dimensional boundary. The gauge 
theory is obtained by assuming that the full SU {N) matrices are decomposed into the 
clusters of submatrices of SU (n) as we have argued. Since the eigenvalue distribution 
is essentially four-dimensional, the resulting low energy effective theory must be a 
four-dimensional field theory. If we further assume A/" = 4 supersymmetry in four 
dimensions which must be present due to the conformal symmetry of AdS^ x S^, 
we may conclude that the low energy effective theory for such a background of IIB 
matrix model is J\f = 4 super Yang- Mills theory. 
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